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Abstract. Using the formalism of general relativity for 
axially symmetric gravitational fields and their sources - 
rotating compact stars - a perturbation theory with re- 
spect to angular velocity is developed and physical quan- 
tities such as mass, shape, momentum of inertia and total 
energy of the star are defined. The change of the inter- 
nal structure of the star due to rotation has been investi- 
gated and the different contributions to the moment of in- 
ertia have been evaluated separately. Numerical solutions 
have been performed using a two-flavor model equation 
of state describing the deconfinement phase transition as 
constrained by the conservation of total baryon number 
and electric charge. During the spin down evolution of 
the rotating neutron star, below critical values of angular 
velocity a quark matter core can appear which might be 
detected as a characteristic signal in the pulsar timing. 
Within the spin-down scenario due to magnetic dipole ra- 
diation it is shown that the deviation of the braking index 
from n — 3 could signal not only the occurrence but also 
the size of a quark core in the pulsar. A new scenario is 
proposed where, due to mass accretion onto the rapidly 
rotating compact star, a spin-down to spin-up transition 
might signal a deconfinement transition in its interior. 

Key words: dense matter - stars: evolution - stars: inte- 
riors - stars: neutron - stars: rotation 



1. Introduction 

In many recent astrophysical applications of the theory 
of dense matter it is necessary to investigate the prop- 
erties of rapidly rotating compact objects within general 
relativity theory. The reason for this development is the 
hope that changes in the internal structure of the dense 
matter, e.g. during phase transitions, could have observ- 
able consequences for the dynamics of the rotational be- 
havior of these objects. Particular examples are the ob- 
servations of glitches and postglitch relaxation in pulsars 
which are discussed as signals for superfluidity in nuclear 
matter (Pines & Ali Alpar 1992) and the suggestion that 
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the braking index is remarkably enhanced when a quark 
matter core occurs in the centre of a pulsar during its 
spin-down evolution (Glendenning et al. 1997). Further 
constraints for the nuclear equation of state come from 
the observation of quasi-periodic brightness oscillations 
(QPO's) in low-mass X-ray binaries which entail mass 
and radius limits for rapidly rotating neutron stars, see 
|Lamb et al. (1998)| , |Li et al. (1999)[ 



The problem of rotation in the general relativity the- 
ory was and remains one of the central and complicated 
problems (Glendenning 1997). Compared to the modern 
methods of numerical solutions to this problem (Friedman 
et al. 1986, Salgado et al. 1994) the method of perturba- 
tion theory (Hartle 1967, Sedrakian & Chubarian 1968) is 
a physically transparent and systematic approach to the 
solution of the problem of stationary gravitational fields 
and their sources. This approach has been applied success- 
fully in general relativity as well as in alternative theories 
of gravitation (Grigorian & Chubarian 1985). 

From a practical point of view in the definition of the 
integral characteristics of the astrophysical objects, it is 
important to analyze the asymptotical expansion of the 
metric tensor at large distances from the stars, to be able 
to compare the results with observational data. One can 
of course introduce the physical parameters of the config- 
uration using the symmetry properties of the object and 
the gravitational field by expressing them in terms of con- 
served quantities. In this work we use the definition of 
the J z projection of the angular momentum of the non- 
spherical rotating star (z is the axis of the star's rotation 
and symmetry of the gravitational field) as a conserved 
integral of motion. It is a well known integral of the non 
diagonal element of the energy-momentum tensor in the 
frame of spherical coordinates. 

Using the method of perturbation theory, we are go- 
ing to calculate the total mass, angular momentum and 
shape deformation from the solution of the gravitational 
field equations in case of hydrodynamical, thermodynami- 
cal and chemical equilibrium for a given total baryon num- 
ber and angular velocity fl of the object. According to 
[Bcdrakian fc Chubarian (1968) the expansion parameter 



is of the order of the ratio of the rotational and gravita- 
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tional energy. It has been shown that, for compact objects 
in the stationary rotating regime without matter flux, the 
first two terms of the series expansion give a sufficiently 
good approximation (see also Glendenning 1997, Weber 
1999). 

The evolution of the rotating stars could have different 
origins and scenarii. Our aim in this work is to discuss pos- 
sible signals for a deconfinement phase transition during 
the evolution of a rotating compact object on the basis of 
solutions for the f2 dependence of the moment of inertia. 

2. Self-consistent set of field equations for 
stationary rotating stars 

2.1. Einstein equations for axial symmetry 

The general form of the metric for an axial symmetric 
space-time manifold is 

ds 2 = e v dt 2 - e x dr 2 - rV[d<9 2 + sin 2 9(d(f> + ujdt) 2 ] , (1) 

written in a spherical symmetric coordinate system in or- 
der to obtain the Schwarzschild solution as a limiting case. 
This line element is time-translational and axial-rotational 
invariant; all metric functions are dependent on the co- 
ordinate distance from the coordinate center r and az- 
imuthal angle 9 between the radius vector and the axis of 
symmetry. 

Reversal symmetry of the time and polar angle <f> re- 
quire that all metric coefficients except to must be even 
functions of the angular velocity 



d4> 

~dt 



(2) 



of the star, the gravitational field of which is described by 
the Eq. ([l]). The physical characteristics of the rotating 
object depend on the centrifugal forces in the local inertial 
frame of the observer. In general relativity, due to the 
Lense-Thirring law, rotational effects are described by ui 
the difference of the frame dragging frequency —uj and the 
angular velocity f2 



uj = + uj{r, i 



(3) 



The energy momentum tensor of stellar matter can be 
approximated by the expression of the energy momentum 
tensor of an ideal fluid 



T l = {e + p)u IJ ,u v -p5£, 



(4) 



where u M is the 4- velocity of matter, p the pressure and e 
the energy density. 

We assume that the star due to high viscosity (ignoring 
the super-fluid component of the matter) rotates station- 
arily as a solid body with an angular velocity f2 that is 
independent of the spatial coordinates. The time scales 
for changes in the angular velocity which we will consider 
in our applications are well separated from the relaxation 



times at which hydrodynamical equilibrium is established, 
so that the assumption of a rigid rotator model is justified. 

Therefore there are only two non-vanishing compo- 
nents of the velocity 



l/ve" — r 2 e^uj 2 sin 2 



(5) 



The equation of state which we will use for our inves- 
tigation of the deconfinement transition in rotating com- 
pact stars will be introduced in Sect. |[ Once the energy- 
momentum tensor (^) is fixed by the choice of the equation 
of state for stellar matter, the unknown metric functions 
v, A, /i, uj can be determined by the set of Einstein field 
equations for which we use the following four combina- 
tions. 

There are three Einstein equations for the determina- 
tion of the diagonal elements of the metric tensor 

G r r — G\ — 8ttG(T^, — T t *) , 
G e e + G* = ^G{T e e +Tl) , 

G r e = 0, (6) 

and one for the determination of the non diagonal element 

G\ = %-irGTl . (7) 

Here G is the gravitational constant and the Einstein 
tensor. 

We use also one equation for the hydrodynamical equi- 
librium (Euler equation) 



H(r, 



= ~ln[u*(r,0)] + const, 
p' + e' 2 



(8) 



where the gravitational enthalpy H thus introduced is a 
function of the energy and/or pressure distribution. 

The parameters of the theory are the angular velocity 
of the rotation and the central energy density e(0) of 
the star configuration. 

2.2. Perturbative approach to the solution 

The problem of the rotation can be solved iteratively by 
using a perturbation expansion of the metric tensor and 
the physical quantities in a Taylor series with respect to 
the angular velocity. As a small parameter for this expan- 
sion we use a dimensionless quantity which is the ratio 
of the rotational energy to the gravitational one for a ho- 
mogeneous Newtonian star E TOt /E glav — (Sl/f2 ) 2 , where 
f^ 2 , = 47rGp(0) with the mass density p(0) at the center of 
the star. This expansion gives sufficiently correct solutions 
already at 0(f2 2 ), since the expansion parameter is lim- 
ited to values O/Oo -€i 1 by the condition of mechanical 
stability of the rigid rotation. This can easily be seen by 
considering as an upper limit for attainable angular ve- 
locities the so called Kepler one n K = yjGM/R% with M 
being the total mass and R e the equatorial radius. For ho- 
mogeneous Newtonian spherical stars fl < = Oq/ V3 
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which is fulfilled in particular also for the configurations 
with a deconfinement transition which we are going to 
discuss later in this paper, see Fig. ||. 

The expansion of the metric tensor is given by 



n ' ' 



(9) 



According to the symmetries of the metric coefficients in- 
troduced in Eq. ([[]) we have even orders j = 0,2, .... for 
the diagonal elements]] 



^l + (Q/n ) 2 ^r,9)]+O(Q 4 ) , 

[i + (n/n ) 2 u(r,6)] + o(n 4 ) , (10) 

and odd orders only for the frame dragging frequency w 

n 



e u(r,d-,n) = e v<-°\r)\ 



e M(r,e ; o) _ r 2r 



w(r,fl ; n) 



O 



q(r,9) + 0(Q 6 ). 



(11) 



The distributions of pressure, energy density and "en- 
thalpy" introduced in Eq. (||) are also included in the 
scheme of this perturbation expansion 

p(r,0;n) = p {0) (r) 

e(r,e-,a) 

H(r,9;Q) 



,(0) 



(r) 



(n/n ) 2 P M(r,e) + o(n 4 ) , 
h(n/n ) 2 H^(r,e) + o(n i ). 



(12) 



All functions with superscript (0) denote the solution of 
the static configuration and therefore they are only func- 
tions of the distance from the center r, the others are the 
corrections corresponding to the rotation. 

This series expansion allows one to transform the Ein- 
stein equations into a coupled set of equations for the 
coefficient functions which can be solved by recursion. 
At zeroth order we recover the nonlinear problem of the 
static spherically symmetric star configuration (Tolman- 
Op penheimer-Volkoff equations), see the next subsection 
2.3 . The first recursion step is to solve Eq. (0) in order 
to obtain the dragging frequency in 0(Q) and to define 
the moment of inertia fo r th e spherically symmetric con- 
figuration. In subsection 2Jj we will consider the second 
order contribution in the £1- expansion (|l^) , ( |l2|) where the 
0(fl 2 ) corrections to the moment of inertia can be found. 
The next terms in the expansion which are of 0(f2 3 ) corre- 
spond to corrections of the frame dragging frequency and 
will be neglected since they go beyond the approximation 
scheme adopted in the present paper. 

2.3. Zeroth order: Static spherically symmetric star 
models 

The functions of the spherically symmetric solution in Eqs. 
( |l(i| ) and (|l2|) can be found from Eq. (||) and Eq. (||) in 
zeroth order of the Q- expansion. 



Notation 



corresponds 



to 



the 



works 



of 



Sedrakian & Chubarian (1968) 



They obey the following equation (Tolman- 
Oppenheimer-Volkoff) 



-G[p(°)(r)+£(°)(r)] 



m(r) + 4np^(r)r 3 



dr ' y " r[r — 2Gm(r)\ 

where m(r) is the distribution of accumulated mass 

m(r) =4tt f e m (r')r' 2 dr' 
Jo 



,(13) 



(14) 



within a sphere of radius r. For the gravitational potentials 
we have 



A (0) (r) = -ln[l-2Gm(r)/r 
i/°>(r) - -\W{Rq)-2G 



r'[r' - 2Gm(r')] 



(15) 
dr' . 
(16) 



i?o is the spherical radius of the star, which is defined 
by P {0) {R ) = 0. The set of Eq.pl) and Eq.(O) fulfils 



the following conditions at the center of the configuration: 
e(°)(0) = e(0) and m(0) = 0. The central energy density 
e(0) is the parameter of the spherical configuration. The 
total mass of the spherically distributed matter in the self- 
consistent gravitational field is Mo(e' o '(0)) = m(Ro). 

2.4- Moment of inertia 

In the first order of the approximation we are solving 
Eq. ( [?]), where the unknown function q(r,6) defined by 
Eq. ( |11| ) is independent of the angular velocity. Using the 
static solutions Eqs.(|l^-(|l5|), and the representation of 
q(r, 9) by the series of the Legendre polynomials, 



q( r ,Q) = ^2 Qm{r) 



m=0 



dP m+1 (cos6) 
dcos9 



(17) 



we find the equations for the coefficients q m (r). It is 
proved that g(r, 9) is a function of the distance r only, 
i.e. that q m (r) = for m > 
[Sedrakian fc Chubarian (1968) 



see 



Hartlc (1967) 



Let us write down the equations for Co{r) — f2(l + 
qo(r)/Qo), which is more suitable for the solution of the 
resulting equation in first order 



l__d 

r 4 dr 



du>(r) 
dr 



4dj{r) 
r dr 



u)(r) = 



(18) 



which corresponds to Ref. Hartle (1967), where it was ob- 
tained using a different representation of the metric. Here, 
we use the notation j(r) = e~( I/(0) ( r ) +A<0) ( r ))/ 2 , for which 
outside of configuration j(r) = 1 , r > Rq. 

By definition, the angular momentum of the star in the 
case of stationary rotation is a conserved quantity and can 
be expressed in invariant form 



J = IT, 



-gdV 



(19) 
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where yJ—gdV is the invariant volume and g — det 
For the case of slow rotation where the shape deforma- 
tion of the rotating star can be neglected and using the 
definition of the moment of inertia Io(r) — Jq(t)/£1 accu- 
mulated in the sphere with radius r, we obtain from Eq. 

,(o) (r)+x (o) (r))/2 u)(r) 



since the expansion coefficients of the deformation A; (r) 
can be calculated from the pressure corrections 



Mr) 



dl {r) _ 8?r 



dr 



= fV[ e C0)( r ) +p (0)( r )] e (- 



(20) 

Using this equation one can reduce the second order 
differential equation (18) to the first order one 

Mr) = 6GJ0M _ (21) 
dr r 4 j(r) 

and solve as a coupled set of first order differential 
equations, one for the moment of inertia ( p0| ) and the other 
( pl| ) for the frame dragging frequency Cu(r). 

This system of equations is valid inside and outside 
the matter distribution. In the center of the configuration 
Iq(0) = and u>(0) = H)q. The finite value wo has to be 
defined such that the dragging frequency u)(r) smoothly 
joins the outer solution 

0(r) = O (l - *^°) . (22) 

at r = i?o, and approaches f2 in the limit r — > 00. In the 
external solution (22) the constant Iq = Io(Ro) is the total 
moment of inertia of the slowly rotating star and Jo = Iq£1 
is the corresponding angular momentum. In this order of 
approximation, Iq is a function of the central energy den- 
sity or the total baryon number only. Explicit dependences 
of the moment of inertia on the angular velocity occur in 
the second order of approximation. 

2. 5. Second order corrections to the moment of inertia 

Due to the rotation in ^-approximation the shape of the 
star is an ellipsoid, and each of the equal-pressure (isobar) 
surfaces in the star is an ellipsoid as well. All diagonal ele- 
ments of the metric and energy-momentum tensors could 
be represented as a series expansion in Legendre polyno- 
mials 



(23) 



g$(r,6) = ^(g fil ,) l (r)Pl(co86) . 
1=0 

It has been shown that the only non vanishing solutions 
obeying the continuity conditions on the surface are those 
with I = 0, 2. 

The deformation of the isobaric surfaces due to the 
rotation can be parametrized by the shift R(r, 9) — r = 
A(r, 9) which describes the deviation from the spheric dis- 
tribution as a function of the radius r in the given polar 
angle 9 and is completely determined by 

R{r,6)=r+(^-\ [A (r) + A 2 (r)P 2 (cos 0)], (24) 



Pi(r) 



dp(°){r)/dr 



(25) 



/ € {0, 2} is the polynomial index in the angular expan- 
sion. The function R(Ro, 9) is the distance of the star sur- 
face from the center of the configuration in the direction 
with the angle 9 to the polar axis. In particular, we can 
define the equatorial radius R e = R(R$ : 9 = tt/2) and 
the polar radius R p = R(Rq ,9 = 0) and the eccentricity 

e=y/l- (R P /Re) 2 . 

The correction to the momentum of inertia AI(r) = 
I(r) — Iq(t) can be represented in the form 



A/ = A/Rodist. + A/shapc + A/Fi. 



Id 



Ah 



Rotation 



(26) 



The first three contributions can be expressed by inte- 
grals of the angular averaged modifications of the matter 
distribution, the shape of the configuration and the grav- 
itational fields, in the form 



AL, 







dh(r)[W { a \r) 



W^ a) (r)/5} 



(27) 



where 

r(Field) 



IE, 



(Shape) 



(r) 



W, 



(Redist.; 



('•) 



n 
n 



(r) = ( — ) {2^,(r)-[/ I (r) + $ I (r)]/2},(28) 

(29) 
(30) 



d Mr) 
dr 

Pi(r) 4 



■ei(r) 



p(o)( r ) +e (o)( r ) 



respectively, which have to be determined from the Eq. 
in second order approximation. The contribution of the 
change of the rotational energy to the moment of inertia 



A/Rotation — ~ 

5 



I0CR0) 



dl (r) [r 2 w 2 (r)e-^« 



(31) 



includes the frame dragging contribution. In this expan- 
sion we neglect the influence of the change of the frame 
dragging frequency, since it corresponds to the next or- 
der of the perturbative expansion ~ 0(£! 3 ). For a more 
detailed description of the solutions of the field equations 
in the ~ 0(fl 2 ) approximation we refer to the works of 
plartle (1967)| as well as Sedrakian fc Chubarian (1968)| . 



3. Model EOS with deconfinement phase 
transition 

For the investigation of the deconfinement phase transi- 
tion expected to occur in neutron star matter at densi- 
ties above the nuclear saturation density no = 0.16 fm~ 3 
several approaches to quark c onfinement dynamics 
have 



been discussed, 



see e.g. 


Blaschke ct al. (1990) 


Drago et al 


(1996) which lead to 



E. Chubarian et al.: Deconfinement transition in rotating compact stars 



5 



interesting conclusions for the properties of quark matter 
at high densities. Most of the approaches to quark decon- 
finement in neutron star matter, however, use a thermo- 
dynamical bag-model for the quark matter and employ a 
standard two-phase description of the equation of state 
(EOS) where the hadronic phase and the quark matter 
phase are modeled separately and the resulting EOS is 
obtained by imposing Gibbs' conditions for phase equilib- 
rium with the constraint that baryon number as well as 
electric charge of the system are conserved (Glendenning 
1992, 1997). Since the focus of our work is the elucidation 
of qualitative features of signals for a possible deconfine- 
ment transition in the pulsar timing, we will consider here 
such a rather standard, phenomenological model for an 
EOS with deconfinement transition. 

The total pressure p({/i;},T) as a thermodynamical 
potential is given by 



where 



(i-x)ph(M,t) 

+XPQ({Hq},T) + p L {{m},T) 



(32) 



h—n,p 



i 



l 



(33) 



is the EOS of the relativistic a — cu mean-field model 
(Walecka model) for nuclear matter (Walecka 1974, Ka- 
pusta 1989), where the masses and chemical potentials 
have to be renormalized by the mean-values of the a— 
and uj— fields m* h = rrih — g^cr, p* h — fih — fjuj&o- The pres- 
sure for two-flavor quark matter within a bag model EOS 
is given by 



P Q (M,T)= J2 P l q d {^T-m q ) - B 

q—u,d 



(34) 



where B denotes the phenomenological bag pressure that 
enforces quark confinement and the transition to nuclear 
matter at low densities. For our numerical analyses in the 
present work, we assume here a value of B = 75 MeVfm~ 3 
which allows us, e.g., to discuss a neutron star of 1.4 solar 
masses with an extended quark matter core. 

In a neutron star, these phases of strongly interacting 
matter are in f3— equilibrium with electrons and muons 
which contribute to the pressure balance with 



p £ ({«},T)= ]T p\ d ^ h T- >mi ) 



(35) 



l=e~ 



In the above expressions p\ d (fj,i,T; rrii) = T; m,) + 

p~l([j,i, T; mi) is the partial pressure of the Fermion species 
i as a sum of particle and antiparticle contributions de- 
fined by 

pt( f M i ,T;m i ) = Ji Tlnll + cxp^fc))] , (36) 

Jo \ llT ) 



where x ± (k) = yjk 2 + mj/T ± n/T. All the other ther- 
modynamic quantities of interest can be derived from the 
pressure (^2|) as, e.g., the partial densities of the species 



nj{{ni\,T) = 



dp({lM},T) 



(37) 



The chemical equilibrium due to the direct and inverse 
(3- decay processes imposes additional constraints on the 
values of the chemical potentials of leptonic and bary- 
onic species (Glendenning 1997, Sahakian 1995). Only two 
independent chemical potentials remain according to the 
corresponding two conserved charges of the system, the 
total baryon number Nb as well as electrical charge Q 



V 



(l-x)n H (M,T) + xn Q ({ N },T) (38) 



o = ^ = (i 



x)QH({^h}, T) + X<lQ({fJ-q}, T) 
q L ({m},T). (39) 



The deconfinement transition is obtained following the 
construction which obeys the global conservation laws and 
allows one to find the volume fraction of the quark matter 
phase x — Vq/V in. the mixed phase where Ph({hii}, T) = 
PqdfJ-q}, T), so that at given ns and T the total pressure 
under the conditions fl38|), ( |3S| ) is a maximum, see Glen- 
denning (1992, 1997). 




Fig. 1. Model equation of state for the pressure of hy- 
brid star matter in /3-equilibrium as a function of the 
baryon number density. The hadronic equation of state is 
a relativistic mean-field model (a — u>), the quark mat- 
ter one is a two-flavor bag model with B — 75 MeV 
fm~ 3 . For comparison the relativistic mean-field EoS of 
|Glendcnning (1989) including p mesons, hyperons and 
muons (incompressibility K = 300 MeV, dotted line) 
and that of Glendenning (1992) including a deconfinement 
transition to three-flavor quark matter in a bag model with 
B = 180 MeV 4 (dashed line) are shown. 
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In Fig. [I] we show the model EoS with deconfinement 
transition as described above. Note that in the density re- 
gion of the phase transition there is a monotonous increase 
of the pressure which gives rise to an extended mixed 
phase region in the compact star after solution of the equa- 
tions of hydrodynamic stability (jl3|). For comparison, the 



relativistic mean- field EoS of Glcndcnning (1989) includ- 
ing g mesons, hyperons and muons (incompressibility K = 

in 



300 MeV, dotted line) and that of |Glendcnning (1992) 
eluding a deconfinement transition to three-flavor quark 
matter in a bag model with B = 180 MeV 4 are sh own, 
see also the monographs by |Glendcnning (1997) and 



Weber (1999 
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Fig. 2. Composition of hybrid star matter in 0- 
equilibrium as a function of baryon number density for 
the model EoS used in this work. 

In Fig. U we show the composition of the hybrid star 
matter as a function of the total baryon density at T = 0. 
Solving the Tolman-Oppenheimer-Volkoff equations (|l3|)- 
( |l5| ) for the hydrodynamical equilibrium of static spher- 
ically symmetric relativistic stars with the above defined 
EOS, we find that a configuration at the stability limit 
could have a quark matter core with a radius as large as 
~ 75% of the stars radius. 

What implications this phase transition for rotating 
star configurations might have will be investigated in the 
next section by applying the method developed in Sect. ^| 
for the above EOS. 

4. Results and discussion 

The results for the stability of rotating neutron star con- 
figurations with possible deconfinement phase transition 
according to the EOS described in the previous section are 
shown in Fig. ||, where the total baryon number, the total 
mass, and the moment of inertia are given as functions 



of the equatorial radius (left panels) and in dependence 
on the central baryon number density (right panels) for 
static stars (solid lines) as well as for stars rotating with 
the maximum angular velocity O max (dashed lines). 




5 

n(0)/n. 



Fig. 3. Baryon number N, mass M, and moment of inertia 
J as a function of the equatorial radius (left panels) and 
the central density (right panels) for neutron star config- 
urations with a deconfinement phase transition. The solid 
curves correspond to static configurations, the dashed ones 
to those with maximum angular velocity fix- The lines be- 
tween both extremal cases connect configurations with the 
same total baryon number Nb/Nq — 1.3, 1.55, 1.8, 2.14. 

The dotted lines connect configurations with fixed to- 
tal baryon numbers Nb/Nq = 1.3,1.55,1.8,2.14 and it 
becomes apparent that the rotating configurations are less 
compact than the static ones. They have larger masses, 
radii and momenta of inertia at less central density such 
that for suitably chosen configurations a deconfinement 
transition in the interior can occur upon spin-down. 

In order to demonstrate the consistency of our pertur- 
bative approach, we show in Fig. |] the values of the ex- 
pansion parameter for the maximally attainable rotation 
frequencies of stationary rotating objects (f2#-/f2 ) as a 
function of the central density characterizing the configu- 
ration. 

In Fig. |5| we show the critical regions of the phase tran- 
sition in the inner structure of the star configuration as 
well as the equatorial and polar radii in the plane of angu- 
lar velocity f2 versus distance from the center of the star. 
It is obvious that with the increase of the angular velocity 
the star is deforming its shape. The maximal eccentricities 
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n(0)/n o 

Fig. 4. Expansion parameter for the maximally attain- 
able rotation frequencies of stationary rotating objects 
(£l/<-/^o) 2 as a function of the central density character- 
izing the configuration. Vertical dashed lines indicate the 
density band for which a mixed phase occurs, see Fig. ^| 
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Fig. 5. Phase structure of rotating hybrid stars in equa- 
torial direction in dependence of the angular velocity f2 
for stars with different total baryon number: Nb/Nq — 
1.3, 1.55, 1.8. Dahsed vertical lines indicate the maximum 
frequency VLk for stationary rotation. 



of the configurations with Nb = 1-3 Nq, Nb — 1-55 Nq 
and N B = 1.8 N Q are e{Q. K ) = 0.7603, e{VL K ) = 0.7655 
and e(H,K) — 0.7659, respectively. Due to the changes of 
the central density the quark core could disappear above 
a critical angular velocity. 

In Fig. |^ we display the dependence of the moment 
of inertia on the angular velocity for configurations with 
the same total baryon number Nb = 1-55 Nq together 
with the different contributions to the total change of 



£ 100 




density 

— shape 

— fields 

rotation energy 



3 

n [kHz] 



Fig. 6. Contributions to the dependence of the moment 
of inertia on the angular velocity. 



the moment of inertia. As it is shown the most impor- 
tant contributions come from the mass redistribution and 
the shape deformation. The relativistic contributions due 
to field and rotational energy are less important. In the 
same Fig. ^ we show the decrease of the spherical moment 
of inertia due to the decrease of the central density for 
high angular velocities which tends to partially compen- 
sate the further increase of the total moment of inertia for 
large f2. There is no dramatic change in the slope of 
at f2 crit = 2.77 kHz. 
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with phase transition 
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Fig. 7. Moment of inertia as a function of angular veloc- 
ity with (solid lines) and without (dashed lines) decon- 
finement phase transition for fixed total baryon number 
Nb/Nq = 1.3 (left panel), N B /N Q = 1.55 (middle panel), 
Nb/Nq = 1.8 (right panel). 
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Fig. [?] shows the dependence of the moment of inertia 
as a function of the angular velocity. It is demonstrated 
that the behavior of 7(0) for a given total number of 
baryons Nb strongly depends on the presence of a pure 
quark matter core in the center of the star. If the core 
does already exist or it does not appear when the an- 
gular velocity increases up to the maximum value f2 max 
then the second order derivative of the moment of inertia 
I(9l) does not change its sign. For the configuration with 
Nb = 1-55 Nq the critical value for the occurrence of the 
sign change is f2 C rit = 2.77 kHz while for Nb = 1.8 Nq it 
is close to fix = 6.6 kHz. 

In order to point out possible observable consequences 
of such a characteristic behaviour of 1(0,, Nb) we con- 
sider two possible scenarios for changes in the pulsar tim- 
ing: (A) dipolc radiation and the resulting dependence of 
the braking index on the angular velocity as suggested 
by Glcndcnning et al. (1997) and (B) mass accretion onto 



rapidly rotating neutron stars. 



4--1- Magnetic dipole radiation 



Due to the energy loss by magnetic dipole radiation plus 
emission of electron-positron wind the star has to spin- 
down and the resulting change of the angular velocity can 
be parametrized by a power law 



n 



(40) 



where K is a constant and n(fl) is the braking index 

n n „ zrn + rn 2 



n(Q) 



<r 2 



21 + rn 



(41) 



where we used the notation I' = 
(dI(Q,NB)/dQ)\N B =const, with the corresponding 
definition of I" (see also Glendenning et al. 1997; 
Grigorian et al. 1999). 

In Fig. H we display the result for the braking index 
n(Cl) for a set of configurations with fixed total baryon 
numbers ranging from N B — 1-55 7V up to N B — 1.9 Nq, 
the region where during the spin-down evolution a quark 
matter core could occur for our model EOS. We observe 
that only for configurations within the interval of total 
baryon numbers 1.4 < Nb/Nq < 1.9 a quark matter core 
occurs during the spin-down as a consequence of the in- 
creasing central density (see also Fig. ||) and the brak- 
ing index shows variations. The critical angular velocity 
^crit(-Ws) for the appearance of a quark matter core can 
be found from the minimum of the braking index Eq. (fi"l|). 
As can be seen from Fig. ||, all configurations with a quark 
matter core have braking indices n(d) < 3 and braking in- 
dices significantly larger than 3 can be considered as pre- 
cursors of the deconfinement transition. The magnitude 
of the jump in n(f2) during the transition to the quark 
core regime is a measure for the size of the quark core. It 
would even be sufficient to observe the maximum of the 
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Fig. 8. Braking index due to dipole radiation from fastly 
rotating isolated pulsars as a function of the angular ve- 
locity. The minima of n(Q) indicate the appearance/ dis- 
appearance of quark matter cores. 



braking index n max in order to infer not only the onset of 
deconfinement (f2 max ) but also the size of the quark core 
to be developed during further spin-down from the max- 
imum deviation Sn = n max — 3 of the braking index. For 
the model EOS we used a significant enhancement of the 
braking index which does only occur for pulsars with peri- 
ods P < 1.5 ms (corresponding to £1 > 4 kHz), which have 
not yet been observed in nature. Thus the signal seems to 
be a weak one for most of the possible candidate pulsars. 
However, this statement is model dependent since, e.g., for 



the model EOS used by Glendenning et al. (1997), which 
includes the strangeness degree of freedom, a more dra- 
matic signal at lower spin frequencies has been reported. 
Therefore, a more complete investigation of the braking 
index for a set of realistic EOS should be performed. 

4-2. Mass accretion 

A higher spin-down rate than in isolated pulsars might be 
possible for rotating neutron stars with mass accretion. In 
that case, at high rotation frequency the angular momen- 
tum transfer from accreting matter and the influence of 
magnetic fields can then be neglected (Shapiro & Teukol- 
sky 1983) so that the evolution of the angular velocity is 
determined by the dependence of the moment of inertia 
on the total mass, i.e. baryon number, 



n 
9. 









AN B ) 



Nb_ 

J=const Nb 



(42) 
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where J = I f2 = const has been assumed. 




Fig. 9. Total baryon number dependence of the spin-down 
rate fi/Q in units of the baryon number accretion rate 
(Nb/Nb) and the corresponding angular velocity (lower 
panel) for different (conserved) total angular momenta 
J[M km 2 kHz] = 100, 200, . . . , 1400. The suggested signal 
for a deconfinement transition in rapidly rotating neutron 
stars with baryon number accretion is a transition from a 
spin-down to a spin-up regime, i.e. a zero in the spin-down 
rate. 



In Fig. [)] we consider the change of the pulsar tim- 
ing due to mass accretion with a constant accretion rate 
Nb /Nb for fixed total angular momentum as a function of 
the total baryon number. Here the change from spin-down 
to spin-up behaviour during the pulsar evolution signals 
the deconfinement transition. When the pulsar has devel- 
oped a quark matter core then the change of the moment 
of inertia due to further mass accretion is negligible and 
has no longer influence on the pulsar timing. However, 
in real systems the transfer of angular momentum from 
the accreting matter can lead to a spin-up already. Then 
the transition to the quark matter core regime should be 
observable as an increase in the spin-up rate. 

It will be interesting to investigate in the future 
whether e.g. low-mass X-ray binaries (LMXBs) with mass 
accretion might be discussed as possible candidates for 
rapidly rotating neutron stars for which consequences of 
the transition to the quark core regime due to mass ac- 
cretion might be detected. Recently, quasi-periodic bright- 
ness oscillations (QPO's) with frequences up to ~ 1200 Hz 
have been observed (Lamb et al. 1998) which entail new 
mass and radius constraints for compact objects. Note in 
this context that the assumption of a dcconfmed matter 
interior of the compact star in some LMXBs as e.g. SAX 



J1808. 4-3658 (Li et al. 1999) seems to be more consistent 
than that of an ordinary hadronic matter interior. 

5. Conclusions 

It has been demonstrated, through the example of the 
deconfinement transition from hadronic to quark matter, 
that the rotational characteristics of neutron stars (brak- 
ing index, spin-down rate) are sensitive to changes of their 
inner structure and can thus be investigated in order to 
detect structural phase transitions. 

The theoretical basis for the present work was a pertur- 
bation method for the solution of the Einstein equations 
for axial symmetry which allows us to calculate the contri- 
bution of different rotational effects to the change of the 
moment of inertia. This quantity can be used as a tool for 
the investigation of the changes in the rotation timing for 
different scenarios of the neutron star evolution. 

The deviation of the braking index from the value n = 
3 (magnetic dipole radiation) as a function of the angular 
velocity has been suggested as a possible signal for the 
deconfinement transition and the occurrence of a quark 
matter core in pulsars. We have reinvestigated this signal 
within our approach and could show that the magnitude 
of this deviation is correlated with the size of the quark 
core, since the influence of the Ae phase crust on such 
processes is negligible. 

For neutron stars with mass accretion, we have sug- 
gested that under the assumption of total angular mo- 
mentum conservation a flip from spin-down to spin-up be- 
haviour signals the appearance of a quark matter core. A 
more detailed investigation is necessary to identify pos- 
sible candidates of rotating compact objects with mass 



accretion (see e.g. Lamb et al. (1998)) for which the sug- 



gested deconfinement signal could be relevant. 

Although the quantitative details of the reported de- 
confinement signals are quite model-dependent and might 
change when one uses more realistic equations of state, e.g. 
including the strangeness degree of freedom (Glendenning 
et al. 1997, Weber 1999), the relation between the magni- 
tude of the effect and the size of the quark core which has 
been found here is expected to be model-independent and 
should be confirmed by subsequent studies. 

Acknowledgments 

We thank A. Drago, F. Weber and N. Glendenning for 
their comments on this work. D.B. acknowledges discus- 
sions during the workshop on "Understanding Deconfine- 
ment in QCD" at the ECT* in Trcnto. Our thanks go to 
R. Avakian, A. Sedrakian and D. Sedrakian for their stim- 
ulating interest in this work and for their useful remarks. 
The work of E.C. and H.G. was supported by the Volk- 
swagen Stiftung under grant No. 1/71 226. G.P. received 
support from the Deutscher Akademischer Austauschdi- 
enst (DAAD) and H.G. acknowledges a stipend from the 



10 



E. Chubarian et al.: Deconfinement transition in rotating compact stars 



Deutsche Forschungsgemeinschaft (DFG), grant no. 436 
ARM 17/1/00. 



References 

Blaschke, D., Kampfer, B., & Towmasjan, T. 1990, Yad. Fiz. 
52, 1059 

Blaschke, D., Grigorian, H., Poghosyan, G., Roberts, C. D., & 

Schmidt, S. 1999, Phys. Lett. B 450, 207 
Drago, A., Tambini, U., & Hjorth- Jensen, M. 1996, Phys. Lett. 

B 380, 13 

Friedman, J. L., Ipser, J. R, & Parker, L. 1986, ApJ 304, 115 
Glendenning, N. K. 1989, Nucl. Phys. A 493, 521 
Glendenning, N. K. 1992, Phys. Rev. D 46, 1274 
Glendenning, N. K., Pei, S., & Weber, F. 1997, Phys. Rev. 
Lett. 79, 1603 

Glendenning, N. K. 1997, Compact Stars, New York: Springer 
Grigorian, H., & Chubarian, E. V. 1985, Afz 23, 177 
Grigorian, H., Hermann, M., & Weber, F. 1999, Part, and Nucl. 
30, 156 

Hartle, J. B. 1967, ApJ 150, 1005; 

Hartle, J. B., & Thorne, K. S. 1968, ApJ 153, 807 
Kapusta, J. I. 1989, Finite temperature field theory, Cam- 
bridge: Cambridge University Press 
Lamb, F. K., Miller, M. C, & Psaltis, D. 1998, Rapid X-Ray 
Variability of Neut ron Stars in Low-Mass Binary Systems, 
and references therein 



ELStro-ph/9802089 



Li, X.-D., Bombaci, I., Dey, M., Dey, J., & van den Heuvel, 

E.P.J. 1999, Phys. Rev. Lett. 83, 3776 
Pines, D., & Ali Alpar, M. 1992, Structure and Evolution of 

Neutron Stars, New York: Addison Wesley, p 7 
Sahakian, G. S. 1995, Physics of Neutron Stars, Dubna: JINR 

press (russ) 

Salgado, M., Bonazzola, S., Gourgoulhon, E., & Haensel, P. 

1994, A&A 291, 155 
Sedrakian, D. M., & Chubarian, E. V. 1968, Afz 4, 239; 551 
Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White 

Dwarfs, and Neutron Stars, New York: Wiley, chap 15 
Walecka, J. D. 1974, Ann. Phys. (NY) 83, 491 
Weber, F., Glendenning, N. K., & Weigel, M. K. 1991, ApJ 

373, 579 

Weber, F. 1999, Pulsars as Astrophysical Laboratories for Nu- 
clear and Particle Physics, Bristol: IoP Publishing 



